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METHOD AND APPARATUS FOR FREQUENCY ESTIMATION USING ITERATIVE 
FILTERING IN A GSM COMMUNICATIONS SYSTEM 

TECHNICAL FIELD 

[0001] The present invention relates to frequency estimation at a mobile station in a GSM 

transmission, and more particularly, using multiple sampled subsets of the frequency control 
channel tone in order to estimate frequency. 

BACKGROUND 

[0002] The most common standard for mobile communications in the world is the Global 

System for Mobile telecommunications (GSM). In one specific implementation, GSM 
utilizes two bands of 25 MHz, which have been set aside for system use. The 890-915 
MHz band is used for subscriber to base station transmissions (reverse link), and the 935- 
960 MHz band is used for base station to subscriber transmissions (forward link). The 
GSM protocol uses frequency division duplexing and time division multiple access (TDMA) 
techniques to provide base stations with simultaneous access to multiple users. 
Transmissions on both the forward and reverse link are made at a channel data rate of 
270.833333 Kbps, using binary Gaussian minimum shift key (GMSK) modulation. 
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[0003] In the GSM protocol, there are traffic channels and control channels. The traffic 

channels carry the digitized voice or user data. One of the control channels is what is 
known as the frequency correction channel (FCCH), which is a special data burst which 
occupies time slot 0 for the very first GSM frame and is repeated every ten frames within a 
control channel multiframe. The FCCH burst allows each mobile station to synchronize its 
internal frequency standard (local oscillator) to the exact frequency of the base station. 

[0004] The data burst carried by the FCCH is nominally at a frequency one-quarter of the 

channel data rate, i.e., 270.833333 -s- 4 or 67.708 KHz. Thus, the frequency correction 
channel is a single tone at the nominal frequency of 67.708 KHz. However, because of 
various factors, such as co-channel interference, fading, and Gaussian noise, the frequency 
of the received FCCH tone may vary from the nominal 67.708 KHz. In order for the mobile 
station to operate optimally, it is important to precisely determine the frequency of the 
FCCH tone to within 1 00 Hz. 

[0005] One prior art method of estimating the frequency of the FCCH tone is disclosed in 

my U.S. Patent No. 5,761,250 entitled "Iterative Filtering Frequency Estimator and 
Estimation Method". In this method, the FCCH burst is iteratively filtered to determine a 
pole estimate. Using the pole estimate, the frequency of the FCCH burst can be estimated. 
Still, the accuracy of the frequency estimation using this technique may not be sufficient. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0006] The foregoing aspects and many of the attendant advantages of this invention will 

become more readily appreciated as the same becomes better understood by reference to 
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the following detailed description, when taken in conjunction with the accompanying 
drawings, wherein: 

[0007] Figure 1 is a schematic diagram of an apparatus for frequency estimation formed in 

accordance with the present invention. 
[0008] Figure 2 is a flow diagram illustrating the method of the present invention. 

[0009] Figure 3 shows a FCCH tone being sampled at the channel data rate in accordance 

with the present invention. 
[0010] Figure 4 shows a FCCH tone that has a frequency offset that results in the sampled 

subsets having a sinusoidal wave. 

DETAILED DESCRIPTION 

[0011] In the detailed description provided below, numerous specific details are provided to 

provide a thorough understanding of embodiments of the invention. One skilled in the 
relevant art will recognize, however, that the invention can be practiced without one or 
more of the specific details, or with other methods, components, materials, etc. In other 
instances, well-known structures, materials, or operations are not shown or described in 
detail to avoid obscuring aspects of the invention. 

[0012] Reference throughout this specification to "one embodiment" or "an embodiment" 

means that a particular feature, structure, or characteristic described in connection with the 
embodiment is included in at least one embodiment of the present invention. Thus, the 
appearances of the phrases "in one embodiment" or "in an embodiment" in various places 
throughout this specification are not necessarily all referring to the same embodiment. 
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Furthermore, the particular features, structures, or characteristics may be combined in any 
suitable manner in one or more embodiments. 
[0013] Turning first to Figure 1, the apparatus of the present invention is shown. The 

frequency estimator 101 comprises an optional bandpass filter 103, a sampler 105, an auto- 
regressive filter 107, a correlator and summer 109, an update calculator 111, and a 
frequency calculator 113. As will be set forth in greater detail below, the data from the 
frequency control channel (FCCH) may be optionally bandpass filtered by the bandpass 
filter 103. 

[0014] The signal x(n) is then provided to a sampler 105 that will form a plurality of subsets 

based upon the signal x(n). The subsets are denoted as Xk(m), where k is the index of the 
subset. Next, the auto-regressive filter 107 filters each of the sampled subsets. The output 
of the auto-regressive filter 107 is then provided to the correlator and summer 109 which 
performs the correlation and summing operation. 

[0015] The update calculator 111 then updates the pole position of the auto-regressive 

filter. The auto-regressive filter 107, the correlator and summer 109, and the update 
calculator 1 1 1 operate iteratively for a predetermined amount of iterations until the pole 
position converges. Finally, a frequency calculator 113 calculates the estimated frequency 
using the output of the correlator and summer 109. 

[0016] The components shown and described in Figure 1 may be discrete processor 

components. However, more likely, the components are more typically implemented as one 
or more digital signal processors that are programmed to perform the indicated function. 

[0017] Next, turning to Figure 2, the present invention will first be described in general 

terms and then the mathematical basis for the present invention will be described after. 
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First, at box 201, the received frequency control channel data may be optionally filtered to 
remove out of band noise and interference. 
[0018] Then, at box 203, the signal x(n) is then sampled into multiple subsets x k (m). The 

number of subsets is designated as p. Each subset Xk(m) contains a series of data samples 
that are taken every p* element in the sequence x(n), starting at index k. Thus, for p=4, the 
subset xi(m) consists of element numbers 1, 5, 9, 13, etc .... of x(n). Further, for p=4, the 
subset x 2 (m) consists of element numbers 2, 6, 10, 14, etc ... . of x(n). The subset x 3 (m) 
consists of element numbers 3, 7, 11, 15, etc .... of x(n). The subset x 2 (m) consists of 
element numbers 4, 8, 12, 16, etc ... . of x(n). In one embodiment, p=4 and there are four 
subsets calculated. In an alternative embodiment, p=8 and there are eight subsets 
calculated. 

[0019] Next, at box 205, each of the subsets are filtered using an auto-regressive one pole 

filter. The pole position of the filter is designated as a(k). After filtering, each of the 
subsets are then correlated and summed at box 207. As further seen below, this results in 
the calculation of the parameter r. Then, the pole position a(k) is calculated based upon the 
arctangent of the parameter r at box 209. 

[0020] The process of boxes 205-209 is then repeated in an iterative fashion N amount of 

times. Each iteration refines the calculation of the auto-regressive filter parameter a(k). 
Finally, at box 211, the carrier frequency is estimated based upon the latest version of the 
parameter r. With the method of the present invention described in general terms, next 
presented is the mathematical specifics with respect to a methodology of the present 
invention. 



[38483-801 4-USOOOO/801 4 App.doc] 



-5- 



11/24/03 



[0021] As noted above, the present invention uses an iterative filtering method after the 

received FCCH signal is first down-sampled into multiple subsets. The FCCH signal is a 
single tone with a frequency unknown to the mobile station (the subscriber) and includes 
156 samples. As seen in Figure 3, the received FCCH burst can be represented as a 
sinusoidal wave 301. The sinusoidal wave 301 is then sampled at the channel data rate of 
270.8333 KHz, i.e., four times the nominal frequency of the FCCH burst. In one 
embodiment, every fourth or eighth sample is then grouped as an element in a subset. Thus, 
there are four or eight subsets, with each subset containing every fourth or eighth sample. 

[0022] As a specific example, assume the received FCCH signal is represented as x(n). 

Optionally, x(n) has been bandpass filtered to eliminate out of band noise and interference. 
Then, x(n) is sampled to form multiple subsets having a nominal frequency of 67.708 KHz. 
If x(n) is at symbol rate of 270.83333 KHz, there will be 4 subsets. If x(n) is double 
sampled, there will be 8 subsets. As used herein, the term sampling also means to divide the 
data sequence x(n) into subsets that contain the data points at equidistant phase (time) 
relationship. This is seen graphically in Figure 3. Note that Figure 3 shows the situation 
where the FCCH tone is at the nominal frequency of 67.708 KHz. Figure 4 shows the 
situation where the FCCH tone has a frequency offset from 67.708 KHz. The multiple 
subsets can be represented mathematically as: 

x k (m) = x(p m + k) 

where p=4 and k=l, 2, 3 or 4 at the normal symbol rate; alternatively, p=8, and k=l, 2, . . 
8 for a double sampling rate. As seen in Figures 3 and 4, which shows a sampling rate of 
p=4, there are four subsets. 
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[0023] Next, each subset is filtered by a one-pole auto-regressive filter with parameter 

"a(k)" which is an auto-regressive parameter. In other words, the filter has a pole at a(k). 
This is represented mathematically as: 

yM = x k {ri)^a{k)y k {n-\) 

where n = 0, 1, 2, . . . , 155 (for a total of 156 FCCH data samples) and a(k) is the auto- 
regressive parameter of the k-th iteration. Note also that higher order filters may also be 
used, such as a two or three pole filter. 
[0024] Next, the filtered output of each subset is correlated and all of the outputs all added 

together. In other embodiments, less than all of the subsets are filtered. This is represented 
mathematically as: 

k q 

where q is the sequence number for the subset. For p=4, each subset x k (m) includes 39 
elements (156 data samples form the FCCH divided by 4). Therefore, q ranges from 1-39. 
[0025] The parameter a(k) is then updated as follows: 

a(k) = fie 3 m 

where /? is a number very close to unity (1) and Zr is the angle of r. In other words: 

\ real{r) ) 

[0026] This iterative process is repeated N times. It has been found that with N = 8, this 

supplies adequate frequency estimation. 
[0027] The estimated base station carrier frequency is then estimated as: 

f = J L.*L 

In m 
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where f s is the sampling frequency, i.e., the data frequency of x(n). The parameter m is the 

interval in the correlation. 
[0028] Using the present invention, it is easier to determine the phase from the signal. The 

present invention estimates the frequency error directly instead of indirectly. 
[0029] While the preferred embodiment of the invention has been illustrated and described, 

it will be appreciated that various changes can be made therein without departing from the 

spirit and scope of the invention. 
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